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GENERALIZED INDICIA! FORCES ON DEFORMING 
RECTANGULAR WINGS IN SUPERSONIC FLIGHT 1 


By Harvard Lomax, Fbanklyn B. Fuller, and Loma Sluder 


SUMMARY 

A method is presented for determining the time-dependent 
(low over a rectangular wing moving with a supersonic forward 
speed and undergoing small vertical distortions expressible as 
polynomials involving spanwise and chordwise distances. The 
solution for the velocity potential is presented in a form analogous 
to that for steady supersonic flow having the familiar “reflected 
area" concept discovered by Evvard. Particular attention is 
paid to indicial-type motions and results are expressed in terms 
of generalized indicial forces. Numerical results for Mach 
numbers equal to 1.1 and 1.2 are given for polynomials of the 
first and fifth degree in the chordwise and spanwise directions, 
respectively, on a wing having an aspect ratio of 4- 

INTRODUCTION 

One of the basic problems arising in the analysis of wing 
flutter boundaries is the calculation of the aerodynamic forces 
on wings undergoing small but arbitrary spanwise and chord- 
wise distortions. Whon the wing aspect ratio is large (actually, 
when the distance between spanwise nodal lines is large), 
these forces are usually estimated by some strip theory in 
which the loading on each spanwise section is approximated 
from that on a two-dimensional wing having the same chord- 
wise distortion. This report is concerned with low-aspect- 
ratio rectangular wings for which tip effects are important 
and the full three-dimensional theory must be used. 

The exact linearized solution for the forces on thin rectan- 
gular wings (limited, however, to the range where effective 
aspect ratio (fiM 2 — 1 A ) i 0 >1) traveling at supersonic 
speeds has been presented by both Gardner (ref. 1) and 
Miles (refs. 2 and 3) in terms of multiple integrals involving 
arbitrary surface undulations. However, the use of such 
solutions in evaluating, numerically say, the forces induced 
by specific wing distortions still presents some difficulties. 
It is the purpose of this report to discuss certain techniques 
that can simplify the labor involved in these calculations and 
to present numerical tables for the forces induced by a class 
of surface deformations, a class general enough to represent 
the first few mode shapes of rectangular plates. 

Mathematically the problem is to find and analyze a solu- 
tion to the four-dimensional wave equation 

<pxxA <pu — — j Vt’t'= o (la) 

Oo 

(whore Oo is the speed of sound, i' is the time, and x,y,z are 
space coordinates) that satisfies the appropriate boundary 


conditions. The particular form of the solution to be 
analyzed differs from those presented by Gardner and Miles 
but its development is based on the method due to Gardner. 

Hadamard (ref. 4) studied a generalized form of equation 
(la) in which the number of dimensions was arbitrary. TTis 
solutions to these generalized equations are fundamentally 
different, depending on whether the total number of dimen- 
sions is odd or even. In fact, the methods Hadamard de- 
veloped apply directly only to equations for which the total 
number of dimensions is odd. Solutions for the even cases 
(such as eq. (la)) axe determined by a "method of descent”; 
that is, the solution for the next higher odd-dimensioned 
equation is found and then reduced by (made independent 
of) one dimension. It is apparent, however, that such a 
technique is in itself by no means unique. Thus, Hadamard 
found the solution to equation (la) by descending from a 
solution to the equation 

Vix J t‘Pn J t<p,r-\-<pn—~i Vt> «/— 0 (lb) 

<zo 

but there are many other partial differential equations and 
groups of partial differential equations governing a five- 
dimensional ( x,y,z£,t ) space all of which satisfy equation 
(la) in a plane f= constant. Gardner discovered a set of 
equations containing equation (la) in a £= constant plane 
which are simpler than equation (la) in that solutions could 
be found and adapted to the boundary conditions for time- 
dependent motion by methods well known to aerodynamicists 
who have studied the flow about wings in steady supersonic 
flight. This is the essential part of Gardner’s contribution 
and it represents the technique upon wjuch the development 
of the solution presented in this report is based. Actually, 
Gardner first applied a Lorentz transformation to equation 
(la) and then used his method outlined above. The appli- 
cation of such a transformation is unnecessary and has the 
disadvantage that the resulting coordinates have lost their 
direct physical significance. We will apply Gardner’s 
method of descent directly to equation (la) and then proceed 
to analyze the solutions so obtained. 

In order to simplify the analysis as much as possible, we 
will limit solutions to the plane of the wing, and, further, 
consider only indicial-type boundary conditions; in other 
words, unsteady motions in which the wing attains instan- 
taneously, at the time zero, a certain spanwise and chord- 
wise distortion which is thereafter fixed. It is well known 
that the transient responses to these indicial motions can be 


i Supersedes NAOA TN 3288 by Harvard Lomar, Frantlyn B. Fuller, and Loma Sluder, 1951. 
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used, in a superposition integral, to obtain responses to many 
other types of unsteady motion; in particular, responses to 
the harmonic oscillations of nonrigid wings. 

Finally, the principal interpretation of the results will be 
made in terms of generalized forces, since these can be used 
directly in either flutter or gust studies, and it will be shown 
that the amount of labor required to calculate such forces is 
reduced by using reciprocity relations derived from the 
general theorems presented in reference 5. 


LIST OF IMPORTANT SYMBOLS 


A 

<h 

& In 

B(p, q) 

q) 

C(x h y x ) 

C L 

Cv a 


O m 


n / 

'-'mg 


m) 

h(x,y,{) 

M 

Ap 

2o 


CD 


2o 

Qr 


aspect ratio 
speed of sound 

amphtude of indicial-downwash distribution 
(See eq. (2a).) 

beta function (See eq. (Bl5a).) 
incomplete beta function (See eq. (Bl5b).) 
influence function for effect of side edge (See 
eq. (A10).) 

lift coefficient, — ~ 

2 d S 

indicial lift coefficient due to angle-of-attack 

2)C , J 

change, without pitching, O l =-5— 

" Oa | a _o 

indicial lift coefficient due to pitching for a wing 


rotating about its leading edge, (7^=-^ 


t-0 


pitching-moment coefficient, positive when trail- 
ing edge tends to sink relative to leading edge, 
moment 
yoSc 

indicial pitching-moment coefficient due to 
angle-of-attack change (without pitching) 

d(7 

measured about the leading edge, C„ a =-g~ 

indicial pitching-moment coefficient due to 
pitching measured about the leading edge for 
a wing rotating about its leading edge, Cm,, 
i>C„ 


£>2 




indicial force coefficient (See eq. 


wing chord 
generalized 
(38)0 * 

generalized indicial force coefficient (See eq. 
(37).) 

distance of wing camber line from z=0 plane 
Mach number 

loading coefficient (pressure on the lower surface 
minus pressure on the upper surface divided 
by free-stream dynamic pressure) 

binominal coefficient, ( n \= — — rr 
’ \mj m! (n— m)l 

c§ 

dimensionless rate of pitching, j=r 

Uo 

free-stream dyn ami c pressure, ~ pUo 

A 

generalized coordinate 

generalized force corresponding to the gener- 
eralized coordinate q T 
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S c 
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tm 
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Uo 

w 
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X3,y 3 ,t3 


Xi.ViA 

Xo 

Xi(v) 

a 

p 
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I 

Po 

<p 


v 


CO 


* 


real part of 
V ( 3 — s i) 2 + (y—yiY 
V(g— a:i) 2 +(y+yx) 2 
V(g— 3 i) ! — Pfy— yi) 2 

wing semispan 
wing area 

area of acoustic plan form 

area of reflected acoustic plan form 

aJJ 

time 

c 

x+Mt 

P 

wing kinetic energy 
wing potential energy 
forward speed of wing 



vertical velocity 

Cartesian coordinates, fixed relative to the fluid 
at infinity 

coordinates with origin on center of wing leading 
edge (See fig. 13.) 

coordinates with origin on center of wing leading 
edge at time zero (See fig. 14.) 
x 
c 

Mx+t 

P 

(Xm—njtJ—1?) 

angle of attack (angle between flight path and 
plane of wing), radians 

I 

wing angle of pitch relative to horizontal, posi- 
tive when trailing edge lies below leading 
edge, radians 

coordinate measuring fifth dimension 
free-stream density 
velocity potential 

portion of velocity potential induced by sources 
in acoustic plan form 

portion of velocity potential induced by presence 
of side edge 

potential function in five-dimensional space 


Subscripts 

A,B,C regions in an x, £ plane (See fig. 7.) 
u upper side of wing, z= 0+ 

1 singularity (e. g., source) position 

I, II,..., VIII regions on wing shown in figure 4 
STATEMENT OF THE PROBLEM 


THE GOVERNING EQUATION 

Assuming a wing’s vertical motion is of such a nature 
that the velocities induced in the fluid are small relative to 
the magnitude of the wing’s steady forward motion, the 
normalized form of equation (la) 

<Pxi-\-<Pn~h<Pu— <ptt—0 (lc) 
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where t=aot', can be used as the governing partial differential 
equation of the flow .field. This equation applies to the 
determination of the velocity potential when the body or wing 
in question moves through the fluid, the axes remaining 
fixed with respect to the still fluid infinitely distant from the 
origin. For convenience we place the wing leading edge 
on the y axis at t = 0 and the side edge on the x axis. The 
wing flies at a constant forward (in the negative x direction) 
speed so at subsequent times the leading edge lies along the 
line x=—Mt, where M is the Mach number, and the side 
edge moves along the x axis as shown in figure 1. 


If the derivatives of h with respect to each of the coordinates 
are small, the two middle terms can be neglected and the 
expression for the boundary condition reduces to 


bh dp 


=wjx,y,i') 

1-0 


We wish to simulate a rectangular wing deformed in- 
dicially by bending in the spanwise and chordwise directions. 
For this purpose, on the portion of the 2=0 plane occupied 
by the wing plan form, the vertical velocity, which determines 
the wing shape according to the previous equation, is assumed 
to have the form 





>o 



Fiqueb 1. — Wing in fixed coordinate system. 


THE BOUNDARY CONDITIONS 

The fluid velocity normal to the surface of a solid moving 
in a frictionless fluid must be zero. If the equation of the 
solid's surface is represented by 

Q(x,y,z,t ')= 0 

this boundary condition can be expressed mathematically, 
in terms of the coordinate system used in equation (lc), as 

bO , bip bO . bip bG b(p bQ - 

bt' bxbx by by bz bz 

Consider a thin surface near the z=0 plane. The equation of 
the camber line of this surface can then be expressed in the 
form 

0(x,y,z,t')=z—h(x,y > t ')= 0 



where c is chord length, a ta is a constant and l and n are 
integers >0. 

The expression (x-\-Mt) 1 is used so that for l ^> 0 the tangent 
to the wing camber line at the leading edge is tangent to the 
flight-path angle of the leading edge. Consider, for example, 
the case Z=l, 7i=0. The downwash 


w v ^(x+Mt) 

c 

represents an infinite class of surface shapes having the form 


Kx,y,i)=^-Kx+M^+f(x,y)] ( 2 ) 

where j{x,y) is an arbitrary function and h is, by definition, 
the distance of the wing’s camber line from the 2=0 plane. 
Since, within the accuracy of linearized theory, the solution 
for the flow about the wing depends only upon the value of 
w u (x,y,t), the loading on all the wings represented by the 
above equation is the same. 

Let us inspect the two special cases 

(i )f(x,y)=—3? 

(ii) f(x,y)=0 

For case (i) 

h(x,y,t)=^(2xt+Mt*) 


and the wing is a flat plate pitching at a uniform rate about 
its leading edge which is following the flight path 


{K)lb 


aio MH 2 

2cU 0 


as shown 3 in figure 2. Hence, at time t the tangent to the 
flight path of the leading edge is 

d(Ji) LB /dt' did' 

-Uo c 

The slope of the leading edge of the plate at the same time is 


and, assuming that thickness and lifting effects 1 can be 
separated linearly, the boundary condition for the camber 
line becomes 



and the two slopes are seen to be equivalent. 


bh.btpbh.bipbh 

dz by by bz 


* The r scale In both figures 2 and 3 Is purposely distorted In order to mate the drawings 
dear. Abe3loassumptionnsedlnsetttngupthebonndary-vahieproblem,bymeansof whloh 
the loading was determined, was that the sur&ce of the wing must remain near the i-O plane. 
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Fiqurb 2. — Flat plate pitching at uniform rate about leading edge. 


For case (ii) 

f,Q=<£$j- o (x+Mtf 

and tie Aving is a plate which obtained a sudden parabolic 
camber at t=0, a shape it maintained thereafter as shown 2 
in figure 3. 
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dh 
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dt 


■0 


ot leading edge 


Fiqurb 3. — Plate with parabolic camber. 


their entirety) only for @A>1, so the opposite edge either has 
no effect or one that can be incorporated by simple super- 
position. 

The wave traces divide the wing area into several regions, 
indicated by the Roman numerals, in each of which the 
analytical formulation for the potential is different. Region 
/ consists of that part of the wing where the effect of neither 
the side edge nor leading edge has yet been felt. -In region 
IT, the side-edge influence is acting (the line y=t is the trace 
of the starting cylindrical wave from the side edge y=0) 
but not the leading edge. Region III is the part within the 
starting cylindrical wave from the leading edge, but outside 
the influence of the side edge. This region, and region V, 
are further subdivided for reasons that will appear later. 
Region TV is a compound region; potential there can bo 
found by adding the potentials for regions II and III and 
subtracting the potential for region I. Region V consists 
of the portion of the wing within the spherical wave origi- 
nating at the wing corner. The flow over the part of the 
wing comprising regions VI and VII has reached a steady 
state relative to a point on the Aving, and the potential there 
is just that for the corresponding parts of a rectangular Aving 
Avith the proper doAvnAvash distribution in steady motion. 
Finally, region VIII is again a composite region, its potential 
being the sum of potentials for regions III and VII loss the 
potential for region VI. 

All the regions just listed, with the exception of region 
V, are actually governed by the three- (total) dimensional 
wave equation and the potential therein could be obtained 
by methods applicable to this simpler equation. However, 
in this report we shall present a unified approach and the 
problem Avill be solved by the same method in all regions. 


The problem is linear, so it Avill bo sufficient to determine 
a solution for arbitrary l and n, and then add results for any 
combination of terms as desired. Thus, ' the complete 
boundary conditions to be studied are 


w B (x,y,f)=~\ 


=®Jn 



(2a) 


over the Aving plan form, and, since the loading is zero over 
the remaining portion of the plane 

=0 off the Aving (2b) 

x-0 

i 

since the loading is given by 

A p_ 4 /bgA 

go UoM\i)t / t «o+ 


d<p 

dt 



Figure 4. — Regions used in the analysis of a rectangular wing in 
supersonic unsteady motion. 


SOLUTION FOR THE POTENTIAL 


BEATEW OF BCBCHHOFF’S FORMULA 


Figure 4 shoAvs the Aving plan form on the surface of which 
the potential is required, together Avith the system of axes; 
also, traces in the 2=0 plane of the wave system set up by 
the indicial motion of the Aring are indicated. The wave 
pattern for only two edges is shoAvn; the flight speed is super- 
sonic so the trailing edge has no effect on the velocities in- 
duced over the Aving surface, and the results are. valid (in 


The solutions developed in the subsequent sections are 
more clearly interpretable if they are compared with certain 
knoAvn results that have already been determined for the 
indicial motion of nonlifting rvings Avith symmetrical thickness 
distributions or lifting surfaces Avith all supersonic edges. 
The purpose of this section is simply to reAueAV briefly some 
of these latter results. 
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As in steady-state wing theory, there is a formula for 
time-dependent flows that relates the velocity potential to a 
distribution of time-dependent- sources and doublets over a 
certain region in the wing plane. This formula is due to 
Kirchhoff, and some of its aerodynamic uses are discussed 
in reference 6. Kirchhoff ’s result is immediately applicable 
in the study of unsteady lifting-surface problems when the 
potential can be represented by sources alone, that is, when 
the upper and lower surfaces of the wing do not interact, as 
is the case in regions I, III, and VI of figure 4. 

Kirchhoff's formula for source distributions can be 
written 


where 


v(z,?/,0 J t)=—~j’J ^dzidy! 
s a 


r, 2 =(x— zi)H-(y — yiY 


( 3 ) 


The brackets on w u indicate that the retarded value is to be 
taken 

[w u ]=w u (.Xi,y u t—r 0 ) 


and S a indicates that the region of integration is the acoustic 
plan form corresponding to the event (x,y,0,t). These con- 
cepts are discussed at length in reference 6. 

As has been pointed out, equation (3) holds for each of the 
regions I, III, and VI, but the area of integration S a differs 
considerably from one of these regions to another. Consider, 
for example, the determination of <p for region III, denoted 
ipnr- Part of the boundary of the acoustic plan form S a 
is found by eliminating T between the equation of the 
leading edge, x 1 —~~MT, and the expression 

(x—xO^ (y—yi) 2 = ( t — T) 2 

which gives the outer boundary, at “time” t, of all the dis- 
turbances that, operating at “time” T, can produce an effect 
at the point (x,y). This boundary is the ellipse 

+(r-?/i) 2 =C (4a) 

where 

Mx+t t _x+Mt 

%m — p * "m p 

If the point (x,y) lies 'within the cylindrical wave from the 
leading edge, that is, — t<Ix<It, the ellipse of equation (4a) 
comprises only part of the acoustic plan form, the remainder 
being bounded by so much of the circle 

(x— Zi) s + (y—yi) 2 =P (4b) 

as lies on the wing at time zero. Figure 5 shows the three 
possible acoustic plan forms for points in region III. The 
limits for the three types are 

(i) t>x>0 

(ii) 0 >x>—t/M 

(iii) —t/M>x>—t 



and these correspond to the subregions III a , lilt,, and III 0 
identified in figure 4. Using equation (3), we can write the 
potential in, say, region III a as 


Villa' 


t—t Jx—^p—Cv—m) 1 To 


where 


l rtr 

J jr— Vf 1 — x j Jx— Jfi)* 7^0 


( 6 ) 




GARDNER'S METHOD OF DESCENT 

Equation (lc) governs a four-dimensional x,y,z,t space. 
Our object, of course, is to find for this equation a solution 
that satisfies the boundary conditions in the z—0 plane as 
specified in equations (2a) and (2b). Obviously, we can 
always construct a space of more dimensions governed in an 
arbitrary way except that it must satisfy equation (lc) in 
an x,y,z,t hyperplane. Then, if a solution in this higher 
dimensional space which satisfies equations (2a) and (2b) 
in the x,y,z,t plane can be found, it represents for £ (the addi- 
tional dimension) equal to some constant the solution to oar 
problem. This characterizes the method of descent. It is 
not obvious, of course, that such a method leads to any 
simplification; but, with a proper choice of the governing 
equation for the new space, such a possibility always exists. 

There are examples where various applications of this 
method have proved to be useful. Hadainard’s use of the 
method, mentioned in the introduction, is classical. A simple 
application of his method is the derivation of the velocity 
potential for a source in a two-dimensional supersonic flow 
field. This potential field (which amounts to a step function, 
the step occurring at the Mach wave) is easy to derive if one 
considers a three-dimensional field with a line of sources 
normal to the free stream and uniform in strength. The 
two-dimensional field mentioned above follows immediately 
by descent. 

In other examples the additional dimension is measured 
with imaginary numbers and the additional law for the ex- 
tended space is the requirement that the functional depend- 
ence on the resulting complex variable shall be analytic. 
The method of descending in the latter case is associated 
with the study of analytic continuation. In particular, 
Riesz’s method (discussed in ref. 7) for solving equation (lc) 
illustrates these concepts. 
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Gardner’s method for solving equation (lc) is to define a 
five-dimensional space in which a potential function i p is 
governed by the equations 

'f'u— 'i'a— ^£e=0 (6a) 

fii—fn—'f'a — 0 (6b) 

and show that solutions to equations (6) in this space are 
general enough to contain general solutions to equation (lc) 
in a plane £= constant. We shall, therefore, proceed by ana- 
lyzing these equations and eventually let £ approach a plane 
in which the boundary conditions of equations (2a) and (2b) 
are satisfied. For convenience, the latter plane is taken to 
be the £=0 plane. 

Since equations (6a) and (6b) are linear, a number of 
possibilities exist for the choice of the dependent variable 
y, z, 0, t). Aside from the more obvious choice 
y, z, 0, t)—ip(x, y, z, i), where <p is the velocity potential 
of equation (lc); for example, one could let yp(x, y, z, 0, £)= 
<p z (x, y, z, t) or again, fo{x, y, z, 0, f)=p(x, y, z, t). These 
various choices amount only to relatively minor differences 
in the detailed technique of the subsequent analysis. If, 
in imposing the boundary conditions of equations (2), one 
is to use only source- type solutions for both equations (6a) 
and (6b), the last choice is sufficient. Therefore, set 

*Kx,y,z,Z,t^=v(x,y,z,t) (7) 

Now differentiate equation (6a) with respect'to'z and set 3 

z=0. 

Defining 


equation (6a) can be expressed in the form 

W tt — Wxx— W {£ =0 (9) 

and the boundary conditions for equation (9) are given 
directly by equations (2). Thus on the wing 


W&x,y,t)=M 


bj/ 

bz 


=known function of y, £ on the wing 

i-0 


Further, the boundary conditions for the original problem 
in (x, y, z, £, t) space require that <p be an odd function with 
respect to z, and continuous across the 2=0 plane except 
over the wing plan form. Thus <p must be zero for 2=0 
except over the wing plan form. The continuation of this 
condition into (x, y, z, £, t ) space then implies, according to 
equation (7), that off the wing 


Hence, both the second partial differential equation and its 
boundary conditions are identical in form to the first set 
given by equations (9) and (10), respectively. Applying 
equation (7) to their dual solution, we obtain the desired 
result 

|j| vK*,y,0,£,i)l ~<p{x,y, 0,i) 

for the potential on a rectangular wing (with /3A>1) in 
supersonic unsteady motion. 


THE GENERAL EXPRESSION FOR THE POTENTIAL 

The method outlined in the preceding section will now 
be applied to obtain integral expressions for the potential 
in any region of the rectangular wing shown in figure 4. 
Consider first equation (9) for TF(£, x, t). This equation is 
the same partial differential equation as that which governs 
supersonic steady flow. Further, tho boundary values in 
the £, x, t space are identical to those representing a thin 
planar wing in a steady supersonic flow. Since the Mach 
number in the steady-flow analog is-V2, the equivalent 
plan form of this wing (shown in fig. 6) is a sweptforward 
wing tip having aU supersonic edges (i. e., the component of 
the free-stream velocity normal to all edges is supersonic). 

Since all edges of the equivalent wing plan form aro 
supersonic, the solution for W can be written immediately 


dPT dp 

d£ £»o bz\ 


l*-o 


'W u (x,y,t)=ai n 



(10a) 


and off the wing 


bW 

bt 


l?~o 


<pt(x,y, 0,f)=0 


(10b) 


Assuming equation (9) to have been solved for the bound- 
ary conditions given by equations (10), we return to the 
second of the set of partial differential equations (6), spe- 
cifically, 


From equation (8), it is seen that the solution to equation 
(9) yields the result 


» It can be shown that the solution satisfies the equation 

is { “ [w«>] [Wm] } 
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in terms of “sources” only, their strength being given by 
equation (10a). Thus, by analogy -with the well-known 
results of supersonic wing theory, we have 

(1I) 

r 

whore r is the area on the wing cut out by the forecone from 
the point (£, x, t). The analytic form of W will differ con- 
siderably in each of the three regions above the equivalent 
wing shown in figure 7. 



Figure 7. — Regions in which analytic form of W(f,x,0 differs. 

The value of W given by equation (11) now becomes a 
boundary condition for the solution of equation (6b). Thus, 
over tho portion of the z= 0 plane for which y>0, £>0, 

the variation of — is now known and for y<0, £>0 
oz ,_o 

tho condition =0 applies. (These conditions are 

«-o 

still not sufficient to determine a unique solution unless 
tho further restriction is imposed that the loading falls to 
zero as the edge y — 0 is approached, i. e., as y-»0-f.) Again 
wo observe that these boundary conditions and the partial 
differential equation (6b) are identical to those studied in 
connection with a stationary planar wing in a supersonic 
stream. As shown in figure 7, solutions from the f, x, £ space 
above the {=0 plane are referred to as W A , W B , and W a , 
depending on the relation between x and £ in a t = constant 
piano. Figure 8 shows the five different boundary-value 
problems formed by the various combinations of W A , W B , 


and W c occurring along constant x lines in the x, £ plane 
and the corresponding regions in figure 4 for which each 
applies. Each of these five problems is directly analogous 
to the boundary-value problem encountered in steady-state 
lifting-surface theory, of a planar, rectangular lifting surface 
in a steady supersonic stream. The "leading edges” of 
these analo gous rectangular plan forms lie along the lines 
&=f, £i=V^— z 2 or depending on the value of x, 



(0 *-/; for regions I,H 



\ y\ 


(Hi) m b ,Y b 



<v) -Ml $x<-t-, SI, YU 


Figure 8. — The five different boundary-value problems in 6, yi, z> 

space. 
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and the "side edge” lies along the line y=0. Hence, by 
means of this steady-flow analog, we can immediately write 
the solution to equation (6b) in the form 


*to.0, M— < 12 > 


where only the area of integration u must he discussed. 

Two possibilities exist for the shape of a. First, if the 
point £, y lies to the right of the dashed lines in figure 8, 
which in the analogous steady-flow problem represent the 
traces of the Mach cones from the leading-edge tips, <t is the 
triangular area shown (for region III a ) in figure 9, part (i). 
If however, £, y lies between this line and the side edge, 
y=0, a is the trapezoidal area shown (for region Y & ) in figure 
9, part (ii). The latter is a well-known result used in steady 
supersonic lifting-surface theory and first developed by 
Eward (ref. 8). The division of the five kinds of problems 
illustrated in figure 8 into the final twelve, represented by 
the regions in figure 4, is brought about by the various 
combinations of W A , W B , and W c that can occur in the area a- 
as the point £, y assumes all necessary values on the wing. 

When has been determined, the potential in the physical 
plane is found by equation (7), or, combining equations (11) 
and (12), 


to » JJ 


If: 


■w u (x t +Mti,y i) dx i dti 

V( f ~ *i) s — & ~ (*— xtf 


(13) 




A detailed analysis of equation (13) for a point x, y, t in 
region V a of figure 4 is given in Appendix A, arid a study of 
this analysis enables one to write the results for all regions 
without difficulty. 

INTERPRETATION OF THE RESULTS 

The results of the rather involved analysis given in Appen- 
dix A can be interpreted in terms of the known solutions for 
simpler boundary conditions. These latter solutions have 
already been reviewed in a previous section in which it was 
shown that the potential on a lifting surface with all super- 
sonic edges can be written in the form 

3 . 

From Appendix A it is found that the potential at a point 
on a rectangular lifting surface can always be expressed as 
the sum of two parts 

<p(x, y, 0, t)=<p m {x, y, 0, f)— (x, y, 0, t) (14) 

where 

J (16a) 

s. 

and 

<p m (x,y, ^ JJ G(x u y,) dxidyi (16b) 

s. 

The value of C(x u y,) is given by equation (A10) in Appendix 
A and the areas of integration, S a and S e , are illustrated for 
the various regions I through 'Viil in figure 16. 

Let us first inspect equations (16) in light of their possible 
analogy with the familiar solution for the steady-state, rec- 
tangular lifting surface. If a rectangular wing having ar- 
bitrary twist and camber is placed in a steady supersonic 
flow, the solution for the potential on its surface can also be 
expressed as the sum of two parts 

<p{x, y, 0)=<p m (x, y, 0)—^® (x, y, 0) (16) 

where, if 

Tc= (x — X]) s — P 2 (y—yiY 

<P m (x,y,0)-~ ( 17ft ) 

ft 

and 

^®(x,y,0)=— iJJ (l 7b) 

ft 

These equations can be construed in the following simple way : 
Equation (17a) represents the potential induced at x,y , 0 by a 
distribution of sources over the wing plan form, each source 
having a strength proportional to the local streamwise slope 
of the upper surface. The area Si, as shown in figure 10, is 
the portion of the wing within the Mach forecone from 
x,y, 0. Equation (17b) lias a similar interpretation; it also 
represents a distribution of sources over the wing, oach 
having a strength proportional to the local slope of the upper 
surface. But the area of integration Si is now that portion 
of the wing within the Mach forecone from the point x,—y, 0; 


Figure 9. — Area of integration <r used in equation (12). 
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that is, within the cone which forms a mirror image of the 
physical Mach forecone in the vertical plane containing the 
wing’s side edge. The potential <p m (x,y, 0) represents the 
difference between the potentials for a wing with a vertically 
symmetrical thickness distribution and a surface with no 
thickness having the same shape as the upper surface of the 
nonlifting wing. 

Lot us return now to equations (15). Just as in the steady- 
state case, <p a) (x,y,0,t) represents the potential induced at 
x,y,Q by a distribution of sources (see eq. (3)) over the wing 
plan form, each proportional to the local slope of the wing, 
but now, since the wing is in motion, with the added con- 
dition that they be local slopes at the appropriate time. The 
area S a , shown in figure 11, is just the acoustic plan form 
defined earlier in the discussion of equations (3) and (4). 
Physically, S a represents those points on the wing from which 
disturbances can, at the time t, influence the flow at x,y, 0 . It 
is the generalization, in the stationary coordinate system, of 
the wing area bounded by the Mach forecone. 

The relation between <p a) (x,y,0,t) and <p® (x,y,0,t) is similar 
to that between their steady-state analogs. Thus, again, 
cp & (x,y,Q,t) represents the difference between the potentials 
for an uncambered nonlifting wing and a lifting surface 
having the same shape as the top of the nonlifting wing. A 
more striking s imil arity lies in the relation between S a and S,. 

We have already seen that S a is the acoustic plan form, 
and, as it turns out, S c isthe reflection of the acoustic plan form, 
in the vertical plane containing the side edge (see fig. 12) — a 



Figure 11. — Acoustic plan form for point in region V a in figure 4. 

situation identical to that existing between & and S 2 in the 
steady-state case. (In other words, S a is the acoustic plan 
form for the event x,y,0,t, and S c is the acoustic plan form for 
the event x,—y,0,t.) Physically S 0 represents the portion 
of the wing’s lower surface containing disturbances which can, 
at the time t, influence the flow at x,y, 0 on the wing’s upper 
surface. At this point the similarity between the steady and 
unsteady solutions ends since the influence of the slopes in 



Figure 12. — Reflected acoustic plan form for point in region V a in 

figure 4. 
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the reflected plan form is not the same as it is for the slopes 
in the basic acoustic plan form; the influence in the former 
case now being given by the integral 0(x lt y 1 ) defined in 
equation (AlO). 

One can show, by simply referring the results given in 
equations (15) to a coordinate system fixed on the wing, that 
equations (15a) and (15b) are identical, respectively, to 
equations (17a) and (17b) when they apply to regions VII 
and VI in figure 4; regions in which, for indicial-type motions, 
the flow is steady relative to the wing. Hence, equations 
(15a) and (15b) extend Eward’s “reflected area” concept to 
all parts of a rectangular wing in supersonic unsteady motion. 4 


THE GENERALIZED FORCES 

REVIEW OF LAGRANGE'S EQUATIONS OF MOTION 

In order to define more clearly the subsequent concepts 
and notation, wo will briefly review Lagrange’s equations of 
motion as applied to distorting wings and will examine a 
simple application to a rectangular wing. 

Lagrange’s equations are usually written 


d bP bP bZ7 
di'bq T £>&■ bg> 


Qr] r= 1,2, 


where 

P kinetic energy of the wing 
U potential energy of wing 
Q r a generalized (external) force 
q T a generalized coordinate 


(18) 


In the present application q, is the amplitude at a given time 
of a polynomial measuring h, the vertical displacement of the 
wing’s camber line from the z=0 plane. Thus, relative to an 
x 3 ,y 3 coordinate system that is fixed on the wing, see figure 13 


Kx 3 ,y 3 ,t')=J2 (19) 




The wing’s kinetic energy can be written 

T= JJ | (20) 

s 


where m is the wing mass per unit plan form area. Using 
equation (19), we find 


h JJ Pr(x 3 ,y 3 )P,(x z ,y 3 )m(x3 ) y3)dx s dy 3 '^ 


oq r 


( 21 ) 


The potential energy is usually diffi cult, to evaluate analyt- 
ically. However, it can often be determined experimentally 
(as will be seen) by measuring the frequencies of the free 
vibration modes. For the present assume that the wing is a 
homogeneous plate of constant thickness. The potential 
energy for such a wing can be expressed as (ref. 9) 


JJ {W- 2 (1-m) 


by 3 2 bir, 2 



dxsdy 3 


which leads to the equation 


( 22 ) 


a b 2 P r b 2 P, 


by 3 2 bars 2 " 1 ” 


if 

b 2 P r VP.^ dXidy3 <28, 


1 b 2 P r b 2 P, 

2 bx 3 2 by, 2 " 


bz 3 b?/3 bx 3 b?/3 

where n is Poisson’s ratio, V 2 = b 2 /br 3 2 -f- b 2 /by 3 2 , and 


- n 2(Young’s modulus) (plate thickness) 3 
• 3(1— y 2 ) 

Now, if the. generalized coordinates have been normalized 
so that each measures the amplitude of a free vibration mode, 
all terms in equations (21) and (23) involving the integral of 
the product of P r and P, are zero. Assuming, henceforth, 
such normalization, we can write 


!ZrJJ PI(x 3 ,y 3 ) m (x^y^dx^dyz+D q r IJ{ (TPr) 1 - 

[w wK^DI}**-** 2 ’ • • • 

(24) 


Finally, dividing through by the coefficient of q, and ex- 
pressing a generalized force as the integral over the wing plan 
form of the product of the rth mode shape and the loadings 8 
2 (Ay), induced on the wing by each of the mode shapes 
considered, we find 


q r +q,wl 


i«T,SSPr(x z ,y 3 ) (^f) dxtfhji 

s s \goA 

ffPAx 3 ,y^m(x 3 ,yi)datdy 3 

S 


(26) 


Figure 13. — Wing in moving coordinate system. 


where to, is the frequency of the rth free vibration mode. 


« It Is of further Interest to notice that equation (16b) can be reduced to a double Integral 
Involving tr.(f,N[) by using, for example, the transformations r*=n+Mfi and r—t—h and 
Integrating with respect to r. 


* We win write (Ap),->jo(Ap/go). where jj Is the free-strcom dynamlo pressure. This Is 
possible without a confusion of notation since the generalized coordinates are depressed as 
■ , - and exclude the term fi. 
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If the free-mode frequencies are experimentally deter- 
mined, equations — such as equation (23) — giving the 'wing’s 
potential energy, never have to be evaluated. Further, in 
such cases, equation (25) applies to quite general wing struc- 
tures with varying density. Usually in the application of 
equation (25), one uses the actual frequency av of the free 
mode but, in evaluating the aerodynamic forces, uses an 
analytical expression that only approximates the rth mode 
shape. Let us examine the generalized force term in equa- 
tion (25), taking, for simplicity, only one term of the sum: 

Qr= q o iJJPrfed/s) (— 'jdzidyi (26) 

s 

According to what has gone before, the mode shape poly- 
nomial P, (r 3 ,?/ 3 ) has the form 

F.C*»0-(f)'(?)‘ (27) 

while (Ap/q 0 ), is the loading coefficient corresponding to an 

indicial deflection (see previous section on boundary condi- 
tions) 

'*-**<) » 

which gives a vertical velocity distribution 

Wu=U 0 q,(l) (29) 

Now a generalized indicia! force coefficient can bo defined as 
follows: 

(f)' (yff [(.Ap/qo),Wy 3 (30) 

s 

(The calculation of these quantities will be elaborated 
in the next section.) Since the generalized force Q r is 
intended to apply to any motion, not necessarily indicial, 
it is necessary to apply Duhamel's integral to the indicial 
force coefficient /}J(f'); thus, 

<31) 

As an example, consider now a simple one degree of free- 
dom vibrating plate. The plate is fixed to a wall and 
restrained along its leading edge. The mode shape is as- 
sumed to have the form 

*~*« (?)'(?)* « 

so for a plate with uniform density and thickness 

Equation (25) now becomes 


605 

For this case, we have the generalized indicial force coefficients 

mi, and my, 



Therefore, equation (33) can be written 



THE GENERALIZED INDICIAL FORCE COEFFICIENT 

It is clear from the previous section that a study of the 
dynamic behavior of rectangular wings moving at supersonic 
speeds can be carried out if one can obtain values of the 
generalized force coefficient, as defined by equation 

(30). We will now show how these values can be obtained 
from the solution to the aerodynamic boundary-value 
problem represented by equation (14). 

It was convenient in developing equation (14) to use a 
coordinate system— x,y,z,t — which was fixed in space so that 
the left edge of the wing moved along the x axis as shown in 
figure 1. On the other hand, in studying the dynamic 
problem it was more convenient to use an x 3 ,y s ,z 3 ,t system 
which is fixed on the wing. In order to convert the results 
in one coordinate set to the other, let us first transfer results 
in the x,y,z,t set to the x^y^z^t sot (shown in figure 14) and 
then, finally, transfer to x 3 ,y 3 ,zz,t coordinates. 

* 5**4 +Mt 

*3= A 

*3’ A 



Figure 14. — Transformations from moving to fixed coordinate system. 
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The indicial force coefficient is defined as follows: 





(36) 


In order to transfer the axes from the set shown in figure 1 
to the more convenient set of figure 14, so that mode shapes 
are symmetric or asymmetric about the wing's spanwise 
center line and the force coefficients denoted j)* can be deter- 
mined, we proceed as follows. First, the loading coefficient 
for a wing in the (x,y) system with downwash given by 



(— D*±)(— D" 



is obtained. This loading coefficient can be written as a 
sum: 



Now the quantity /jj is defined in the x^y^z^t system as 


DETAILS OF CALCULATION 

The details of actually evaluating the indicial force co- 
efficients from the solution for the potential presented in the 
first part of this report are discussed in Appendix B. Con- 
siderable labor is involved in such calculations, and an 
attempt was made to discover recursion formulas by means 
of which certain derivatives, for the rectangular wing, could 
be expressed as combinations of others. This attempt was 
successful and yielded the following results: 

Consider equation (36). Integrate the x integral in this 
equation by parts, setting 


u(x)= f y*^—dy; dv(x)= (z+Mi ) hlx 

Jo 2o 

Then, since by equation (B7) in Appendix B 


one finds 


b A y l l A p l Un n 

= : , £^>u 

bx go c g 0 


Fi*=— - 
” j + 1 




(38a) 


This last integral can be written as 

“r 3 P (0 (fr S C) 

By using equation (36) we find 


where all forces are responses to a unit indicial disturbance. 
Note that if equation (37) is applied in the case of a wing 
cantilevered on a wall, both n and g must be even in order 
to satisfy the boundary conditions of reflection in the wall. 

By superimposing boundary conditions and their result- 
ing solutions, one can further show that the value of /}* 
given by equation (37) is valid for all reduced aspect ratios 
/SAL greater than 1 in spite of the fact that the value of F)* 
given by equation (36), as it stands, applies only to wings 
for which PA is greater than 2. 

Given /}£({'), one can determine the generalized force as- 
sociated with the generalized coordinate g r by means of the 
superposition integral as illustrated by equation (34). 


Inspection of equation (37) shows that the same relation 
holds for the generalized indicial force coefficients that is, 

f£=j^i {fo-'-l-ff+M (38b) 

From this relation, it is seen that only the forces FJJ need 
be determined by integration; the forces for higher values of 
the index l can be found by combination of results for dif- 
ferent values of the mode shape index j. 

As a simple illustration of the results presented so far, wo 
can calculate the indicial force derivative for the cases 1= 
n=g= 0, j— 0, 1. The case j=0 corresponds to the indicial 
lift coefficient for a flat, sinking, rectangular wing, and the 
case for j= 1 corresponds to the indicial pitching-moment 
coefficient for the same wing. Since n=g=0, equation (37) 
gives 

f?o=Ko 


Thus, with j=0 and identifying —a w /U 0 as angle of attack 
a, one finds from Appendix B 

_ »{J C0S '‘ M h ±JI T W-W+ 

M+i- u -M = l 

"K'T k) ^12=1 
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Next with j= 1, and using O m J to designate the pitching 
moment measured about the leading edge of the wing, 

j cos-* -w]- 

— K'-m) ^sft 


that is, the lift coefficient for a pitching wing equals the sum 
of the lift and pitching-moment coefficients of a sinking wing 
(primes indicate the wing is pitching about and moments are 
measured about the wing leading edge). Hence, 


M *+ 1 


2 C 1 r/'i . 1 2 \ • Afto — 1 i 

=m{^L( i + 2 , ») C0S 


cos 


3-Mo 




(M- 


-w]}s 


1 <to< 1 


M+l 


M—l 




These expressions agree with those given by Miles in refer- 
ence 2. 

The above results can he used to demonstrate the useful- 
ness of equation (38a). Taking j=n=g= 0, 1=1 in that 
equation gives 

F&=F$-FT 0 


or, for the present case, 


no m> jw) 

Joo— Joo JlO 


which represents the equality 


A further application of equation (38a) provides the pitch- 
ing-moment coefficient for a pitching flat rectangular wing. 
Thus, with l=j= 1, n=g=Q, equation (38a) gives 


which becomes 


and so 


F\°o=m-F%) 


<V 


V M 

2 \ dmfUo 



0 L '—C La +C m 


From equation (B21) in Appendix B it is found that 


Ho Uo 


cos \M— 0%)+ 




Combining, we find 


n t £_ 


■g{ 2+ f‘"’ ^ tS-6Af^.+i^(JV^+3)faU } ; 0<<,<j+j 

2 f lP+lffo’ , 2M , 8— ifto— (JC+2)(,' /m; — ml 

= “H\+— 3 “ ' m -iT'+aV 003 (M “«+ 5 r AM1 _ Mtl) -j- 


i r 
24AL 


24+lM+l' 


-8io — 6(M — 1)U?+(M- 




_u 


1 <to< 1 


M+l 


M—l 




413072— ST 10 
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Another relation among the generalized indicial forces 
can he derived by means of the reciprocity relations, given in 
reference 5. The details of the derivation are given in 
Appendix C and there results 

<-»-(!)£ < 39 > 


Equation (39) can be used in two ways; one, as a means 
for checking the internal consistency of a set of calculated 
generalized indicial forces, and the other, as a means for 
expressing a given force in terms of a set of others. 

Consider, as an example of the former use, the case for 
which l=j= 0. Then 

ft n c=foi 


From equation (37) we can express this relation in terms of 
the calculated quantities thus 




g-bn—u—r 


FZ= 


£ (-D’ 

►-o 






( 40 ) 


The inner sum on the left can be evaluated. Thus one has 


a?=[l-(l-x)p=S (-!>* ( V ) (l-x) 
o 




r<V 

V 


Equating coefficients of x, 
and equation (40) becomes 

CONCLUDING REMARKS 


If now 7i—l, g= 3 the following relation results 

(7^-J1®+| [(A?S-F«)+3(FS-FgO]+ 

3 (jJ(F%-F®+2 (jJ(F^-F^)= 0 

which provides a useful check on the computed quantities. 

Next let us solve equation (39) for a given force. Perform 
the sum operation 

on both sides of equation (39), and reverse the order of 
summation on the left side. There results 


A method is presented for evaluating the generalized forces 
on a rectangular wing flying at supersonic speeds and having 
an aspect ratio such that pA>l. The generalized coordi- 
nates used to define the wing’s behavior are the amplitudes 
of downwash distributions expressed in terms of polynomials 
in x and y, the chordwise and spanwise directions, respec- 
tively. 

Numerical results are presented in table I for generalized 
indicial forces on a -wing having an aspect ratio of 4 and 
flying at a Mach number equal to 1.1 and 1.2; the polynomial 
coverage being 0<Z<1 and 0<7i<5, where w~x l y n . 

Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., June SO, 1954 



APPENDIX A 


EXPRESSIONS FOR THE POTENTIAL 

la order to write the expressions for the potential in all 
regions shown in figure 4, it is sufficient to derive in detail 
only that for region V. Having carried out this analysis, 
one can determine the expressions for potential in other 
regions without difficulty. 

Consider, therefore, equation (13) and let <j and r apply 
to region V a . First, it is necessary to determine the poten- 
tials W A and W B in the space. From equation (11), in 
conjunction with figure 7, it is found that 

n/= _ip+^, ri-vg^qg 
A *J*-VP=i? ‘Jo VG-^-fiM*-*!) 1 

(Al) 


w,=-l f° ^ wJn±MhM 

•a- J X.tt.) 1 J -xjM V(t-i,) 2 — (x-r,)-' 

lf I+Vp "^, ?T Wu { Xl+ Mt u y,)dt x 

^Jo ‘Jo V(i-* i)*-fe*-(a— *0 S 

(A2) 

where 

"With the values of W given in equations (Al) and (A2) 
it is possible now to solve equation (6b) for figure 8 giving 
the required data in the plane. Thus, if R 2 — (£— £0 2 — 

(: V—ViY 


***>— %,r *3_ir 

*VJ+r-< J?+(,-n) ttJ, J ’ r Jf+jf-V< iz i* Jf+(*-n) «- 


9 TFj 


1 r-f+jr+Ve^F* rVe-i* HA— IF, 1 f° IF, 1 P-H-i-jr r< w. 

£ ^ dyA dtA- T f+l\ dyA c/f, %+ 

y j£-(y— yi) */£+(y“yi) ^V/0 */€+(y+yi) 


if *fi=s4ir^* r 

^ J H-y— Vf 3 — ^ */H-(y— yi) ** 0 


IF-TF , 1 f-f-H-Ve-f pVe=F» PF-TP 


H-CH-in) 




A 


(A 3 ) 


Now apply the operation of equation (7) and the potential <pr 0 is given by 



ZiW A 


Rx 3 


-T'd^-f , 

Jy 





pVe-i 1 

- 

Jf-jn 




A, 3 



where A l 2 ={ 1 a — (y—yi) 1 and the bars on the integrals signify 
that the finite part of the integral is to be taken in the sense 
defined 1 in reference 10 and that the order of integration 
cannot, in general, be reversed. 2 

For convenience set 


1 10 

S I* (AS) 

7T 1 

where /„ is the nth integral group on the right-hand side of 
equation (A4). 

Consider the first of these integral sets. Using equation 

1 For tho subsequent analysis to hold, the definition of the finite part given In reference 10 
Is essential. This definition differs from that given by Hadamard when It applies to multiple 
Integrals. 

1 Since the order of Integration plays an Important role In the following development, 
integration first with respect to x and then with respect to y will be denoted fdy fdxflp, y) 
while Integration first with respect to y and then with respect to x will be denoted fdx fdy 
(x, y), When the notation y) dydx Is used, the order of Integration Is immaterial. 


(Al), we can write 

r C J r $i<& f*+Vp=i? 7 

Il ~),-i dyi : t-* dXl 


p-Vcx-r.mc wJjxr+Mt^dt, 

jo ^ V(<— y 2 — (x— Xif — fr 

In order to simplify this expression, the order of these in- 
tegrals will be rearranged so the integration with respect to 
£1 can be carried out first. The technique of changing the 
order of repeated integrals with strong singularities set forth 
in reference 10 will be used here. Consider the change of 
order in the £i,Zi plane. Pretend for the moment, that the 
ti integration has been carried out. Then the highest order 
singularity (since w u is bounded) in the &, Xj plane has the 
order 3/2 which is weak in the sense that no residual occurs 
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wlion the sequence of integration is reversed. The top of 
figure 15 shows the area of integration, so immediately 



[tf-ly-yW 


:-V Ci-iOH-ft 5 


sr 


j{t— #i) s — (x— zi) 2 — fi 4 


To change order in the U plane, consult the bottom of 


figure 15. In this case an inherent singularity exists at the 
confluence of the singularity lines of the integrand; namely, 
where Zi—y—yi and ti=t— ^{x—x) i 2 +£i 2 . The change of 
order can therefore not be performed directly, but account 
must be taken of the existence of a residual term (see ref. 10). 
This residual is defined as the difference between the two 
integrals taken in different orders over a vanishingly small 
region surro unding the inherent singularity (the dotted 
region in bottom of figure 15). The residual is then, 


' fijja 

I^-fcr-ydT 2 , 


rt-ifUT- 

Jt-ro-t 


, w u (x 1 +Mt 1 ,y 1 )dti 

/(f-fOMz-zO-h 2 " 


f wJxx+M^y^dtrf 

Ji-n-t Jr 






where r<?= ( x — ®i) 2 + (y — yi) 2 . The second integral varushes 
(see ref. 10), and, passing to the limit e— >0 in the first integral 
there results 

T> it w u (xi+Ml—Mr 0 ,yd tt [ui„] 

Ki ~ 2 r„ 2 r 0 

where the square brackets again mean that the retarded 
value is to be taken. Thus, the integral I\ can be reduced to 



In the same way, the integral J 2 can be reduced, and 



which is recognized as Kirchhoff's formula, equation (3), 
with an acoustic plan form bounded by the circle 

(z— aj])*+(2/— 2/i)*=* s 


The reduction of the integrals I 3 , /<, I&, and I B is quite 
similar, leading to the sum 



(A7) 


Examination of the limits on these integrals shows their 
total area of integration is that shown in figure 11. But 
this area corresponds exactly to the acoustic plan form S a 
for a point in region V^I Hence, denoting the combination 
of terms in equation (A7) by <p w we can write simply 

isJPs 1 ** (A8) 

It now remains to calculate the integrals / 7 through Jio. 
Designating their total effect on the potential, ip ® ono can 
readily show (since no inherent singularities arise in these 
cases) that 


Figure 15. — Areas of integration used in analysis. 



GENERALIZED INDICIAL FORCES ON DEFORMING WINGS 


611 


,.u> i 



rx-H/p-is+n)* , 
dy i I , dx i 



V4yy t yQdfi 1 p 

[(i— fj) 1 — ro 2 ] ^ Jo 1 J*— 




"'- r ‘ Vi mWufa+MtuyJdt! 1 . p p-r, V4m «^+Mi, yO^i 

o [(t-ti 2 -r<?\ sKt-ti) 3 -^ V s J o 2/1 Jxi(rt*) 1 J^r./Af [(i — «,>* — r 0 ^ 


m 


where r 2 = (x-Xi) 2 -\- (y+y x ) 2 . Now let 

[/: 


t-ri 


■yfiwi w u (x x +Mt u y^dii 


O(xi,y0~-{ 


*iw [(t-ttf-ro 1 ] V(^-y 2 -n 2 


£i<0 


p-o Viy^ ro„fe+Mi, i/i) <fti ^ 

Jo [(f-^-nfl 


(A10) 


In terms of this expression, equation (A9) can he written 
simply 

v>v=^5 jj 0(x l} yi)dxidyi (All) 

Wr a 

where the area (S e )r a is illustrated in figure 12. 

In order to give expressions for the potential in every 
region of the wing shown in figure 4, one can show that it is 
only necessary to vary the areas over which the double 
integration in equations (A8) and (All) are carried out. 
This is evident in connection with tho source portion <p m , 
for in every case 

U (Ai3> 


and only the acoustic plan form S a changes with the region. 
In the case of ip™, the part of the potential due to the exist- 
ence of the side edge of the wing, equation (All) can be 
generalized and written 

p (2) =:3 J I* 0(x i, 1/0 dx x dy x (A13) 

s c 

where the integrands are defined in every case by equation 
(A10) and only the “reflected” acoustic plan form S e 
changes with the region. The region S c is always bounded 
by portions of the “reflected” circle 

(r— Xj) 2 +( y+ yi y=p 

and the “reflected” ellipse 

(m x i ~ Xj *) + ( y+yy =t ™ 

Figure 16 shows sketches of both S 0 and S a for all regions in 
figure 4. The absence of a sketch indicates that the corre- 
sponding integral does not exist for that region. 


APPENDIX B 


THE GENERALIZED INDICIAL FORCES 


THE LOADING COEFFICIENT 


In order to determine total forces acting on the wing, it is 
first necessary to obtain expressions for the loading coefficient 
Ap/q 0 . According to the linear theory 


A p 4 dip 

g 0 ~~ UoM- dt 


(Bl) 


so it is necessary to differentiate each of the expressions for 
potential. As an example, consider, as in Appendix A, just 
region V a of figure 4. The loading coefficient will be divided 
into two parts Ay m jqo and Ap m fq a to correspond to the po- 
tentials <p il) and <p C2) . Thus, using equation (All) 


‘Ap\ ® 4 / f-r+‘ p+Vs-Of+r.)’ ftC 

.q.o) ra~rVoM\J 0 ^J^VFSW dt 


X 

X 


-V+dP-x* 


— JP—x 1 


Jx,( jr+m) or J 


(B2) 


since the derivative passes the x x ,y x integration without 
effect. Referring to equation (A10) for the function C(x x , y x ) 
we next find its derivative with respect to t. Write r=f — t x ; 
then for x x <^0 


v C t+Zi!M yi m w u {x x +Mt~Mr,y x )dT 

’ 1 Jr i (t 2 — To^Vr 2 — r x 


and 


dO_ 
dt 


■yfem Wn(Q ,y) 


('t+xjM^m ^JyW u {x x +Mt~MT,y x ) j- 
Jri (t 2 — To 2 ) V T 2 — Tj 2 ) 


dr (B3) 


Notice that if w u does not depend on (x x -\-Ml x ) the integral 
term in equation (B3) vanishes, while if it does, then the 
integrated term is zero. Next, for £i>0, 


GCr.,*)- ir 

Jn (r— tv) Vt — r x 

and 

ap JSAaJuM ■ ^ 

& (^-ro 2 )#^? Jo 

m 

In this case, both terms exist unless w u is not a function of 
(zx+Afifi), in which case the integral vanishes. 

Substitution of equations (B3) and (B4) into equation 
(B2) will now yield an expression for the loading coefficient 
corresponding to the influence of the side edge; 
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Figure 16. — Sketches of areas of integration, S c and S a , for all regions 

in figure 4. 



Figure 16 — Continued 


/Ap\ 4o„ ■ f ['~ ,+L p+yp-dr-ig.P ^J^ lXl l y 1 n dxi , ,, f~ g+t , p+V.' ! -(r+r ,) 1 ri-n ^ 4 ^ (Xj+MtJ'-'ySdtj 

\ 2 o/ y, ^U 0 Mc ,+n \J 0 ‘Jx-Ve-OrtvO’ (f 2 — ro 2 ) V^~ r i s J° Vx-Vu-Crt-yi ) 1 'Jo [(f— <i) ! — ro^VfZ— -fi) 2 — ri 4 


J -tr+Vu^S f0 


V4yy! xjyj'dx^ 


: [-’"'"-"ij,, r <&, r ” 

Jo Jx-V t»-(y+yi)» Jo [(£ — ii) 2 — ro 2 ] V(^ — ^i) — t’i 


(r+n) 1 (f 2 — ro 2 ) Vf 2 — Ti 

t_ri V4yyi (z 1 +Mfi) ; ~Vi g J£i 


Jo Jx.(K+f.) J-Cn/xo t(t— f]) 2 — r 0 ^ ■sj(t—t\) —rj J 


(B5) 


Tlio explicit form of w u , given by equation (2), bas been 
inserted and it is assumed that Z> 1. 


The portion of the loading coefficient corresponding to 
<p{r l) can be found readily and is 


(?) 


( 1 ) 




In 


'^AfUoc l+tt 


S r' +t u Ax+^-(y-yyy+(x-^-(3J- yi )*y d 
\Jo Vl je-fy-yd* y 


rv+t rx-l 

my yi*dyij^_ 


i+ye-(v-n) ] (£ — ro)] 1-1 


Jp—l F— Fi) 5 


r 0 


dxi — 


AfZ yiJyi , 7 dzi + 

Jo J X — JP - lf - fl ) 1 '0 

mi r ^ v ,-dy, r fa+w-^r- 

Jo J xi(j— in) ro 


(BO) 
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Figure 16 — Concluded 


It is clear that, evon for small values of the indices l and n, 
the required integrations for the determination of total forces 
on the wing pose formidable problems. There is, however, 
a property of the loading coefficient corresponding to vertical 
velocity distributions of the type chosen here (eq. (2)) that 
will materially shorten the requisite labor. This may be 
expressed as follows, adopting the convention that Ap In /g 0 
corresponds to a downwash distribution, proportional to 
(z-\-Mt) l y n : 


or 


d Ap tn l Ap l ~ 1,H 

bx go c go ’ 


Z>0 


(B7) 


f ^-^(x l ,y,f)dx 1 , l > 0 (B8) 

go C J _ir: go 


DETAILS OF EVALUATING THE GENERALIZED INDICIAL FORCES 

In calculating the generalized indicia! forces by means of 
equation (36), it has been shown that only the value zero 
need bo taken for the index l. Thus we must find 


r_: (z+Mtydxjy ^ dy CB9) 

The values of the loading coefficient Ap oa /g 0 are found by 
differentiating the expressions for potential given in the first 
part of this appendix. 

It is convenient, in evaluating equation (B9), to consider 
the integration with respect to y first. Setting 

it is found that L seems to have different representations 
according to the interval in which x lies. These expressions 
can, however, all be expressed by the same formula. The 
portions of L corresponding to the parts ?> <u and <p w of the 
potential are similarly signified, and we have 

(n+g+ 1)! l K ^ n +9)+Kx(n+g)]~ 

2 §(2^) ^+r+i-l } (Bill 

[go (n+g)+K M (n+g)\ ' (Bl2) 

where 


Ko(n+g)=t*+‘+'R.P. 
M 


rqotr-i C-x It) 

.1 sin n+lt+1 OdO 

K M (n^g)=M-t m ^R.P. J^ 1 sin n+rfl OdO 


and [n/2] means the greatest integer contained in n/2. The 
function J(n,g) may be expressed as summations, and it has 
the property 


J(n,g)=J(g,n) 

The sum formula is, with g+p=n 


(B13) 


ig 




P J r2g — 2f4-2j+3 1 


I) 


(B14) 
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Values of the function J (g,n) 

n. n 
9 \ 

0 

1 

2 

3 

4 

5 

0 

x — 2 






1 

i 

1 a. 4 
~4 x+ 3 





2 

5 8 

4 * 3 

1 

2 ' 

29 16 

64 T 15 





11 

21 8 

1 

53 . 32 





6 

64 T ~t5 

3 

256 x '35 




4 

189 32 

64 x 5 

11 

15 

129 128 

256 x 105 

1 

4 

5329 
16384 x 

256 

315 


5 

71 

15 

165 ,64 

256 X ‘ 21 

37 

84 

975 . 64 

4096 x "‘ 63 

1 

5 

11801 512 

65636 x+ 693 


where 



is the binomial coefficient 

(V \ _ Pj 
\2i/ (2i)!(p — 2i)l 


and -B(p, g) is the beta function 

•B(p,g)=J' x)*-’dx 

=2 J* sin 211-1 0 cos 2 * -1 0d0 
,=r(p)r(g)/r(p+ g) 


(Bl5a) 


The function J(g, to) has been calculated for g, to taken 
0, 1, 2, 3, 4, 5. Because of the property (B13), it is only 
necessary to give a triangular array, which appears in the 
above table. 

Now consider the functions Ko(v) and Ku{v), defined after 
equation (B12). It is convenient, for computational pur- 
poses, to express these in terms of the incomplete beta 
functions, defined as 


rcos -1 to 

Bi-Ap/j)=2j a sin 2 *- 1 0 cos 2 * 1 


ede 


-r 


(Bl5b) 


r-'a-Q*" 1 di 


A tabulation of the incomplete beta functions is available in 
reference 11. Note that when the symbol B is written 
without a subscript, the complete integral is meant, that is, 
in equation (Bl5b), x equals 0. It is necessary to exercise 
some care when interpreting Ko(v) and K M {v) as beta 
functions because of the upper limit. Thus, since 


1( X,l) sin’ +1 OdO 


we have the following cases: 
(i) x>t, RJP. cos -1 


(ii) 0<x<t,RJP. cos -1 0=cos _1 j)~ v ~ c 08_1 (j) 

\_ w I)] 

(hi) — t<x<0,RjP. cos -1 ^—0=cos" 1 ^ 

(iv) —Mt<x< — t; RJ?. cos -1 ^ — ^=0 

0 

A simil ar line taken with K M {v) leads to 


Q x>t,K M (v)=0 

® KM- [>W (t 2 - I)] 

SB) 

-<■.».>■ Cir £)] 

Civ) 


The generalized indicia! force can now be expressed as 


Ff g 


&tOn S 1 f J(g,ri) , p, -,s 

TrMU<fi> + * + * + ' 1 4 L 2 r+ * +2C 

['*+”)+*^+")]-§(2”) 


n »lgl ~ | 

(TO+g+l)lJ 

gf+n+1— 2;x 

g+TO+ l—2/i 


(BIG) 
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where 


J c~Mt r f cos -1 (-i/i) 1 

(B17) 

*&(„)« J^" 1 sin’+^^J 

(B18) 

It is convenient to express these forces in terms of dimen- 
sionless quantities. Thus setting 


x t 

Xo=-> <o=- 

c c 


we have 


*J,!(v)=^+'+ s J_ 1 ^'°(x 0 +Mfo) / dx o [V+ 1 12. P. 


^/h(v)=P+'+ 2 


and 


rcos-H-x<jtt) ~l 

J sin' +1 I 0 } (v) (B19) 

cw4f(^r *?. 

rcos -M2±!! *1 

J X,+Mh sin r+1 Odti\=c i+ '+ 2 IJ (v) (B20) 



Fiqtjbb 17. — Areas of integration used in analysis. 

(i) Xo=*o; 0<to 

(ii) Xo—l—MW, <to< ^£—1 
Ciii)Xo=~fo;^ = j<<o 


'S- JSrfe [w+t- 1 ’’ 2 5^1),] [« &+»)+ 

/A\ 0+n ~ il1 

* <*-*}-§ (£) <*-«>+* e*- 1 )]} 

CB21) 

The integrals /^(^) and I{t(v) can be simplified by reversing 
the order of integration. This can be accomplished in a 
straight-forward manner by merely inspecting the region of 
integration in the x a ,0 plane. Consider first the integral 
Jl(v). Depending upon the relation between the chord 
length and the time, we see — from figure 17 — that reversing 
the order of integration results in three different possibilities 
for the upper limit of the 0 integral. However, if we define 
Xo such that 


then, in every case, I$(y) can be written 


Ii(y> 


t 0 r+i rco^H-Xtfo) 


3+1 

U )+r+3 ■ 


I 


skr H ede— 


f+r+3 H-l A’-Ln Fear 1 (-Xfl/W 

2 (-1 y Mi+1 ~ r J 0 sin' +I 6 cos T QdQ 

(B22) 

and, similarly, it can be shown that 

.1+MxjU 


™sin-^+ 


MtZ +,+i 

j+v+2 


it-, 


wr 1 ( — Xo/io) 


sin* 4-1 0 cos r Odd 
(B23) 



APPENDIX C 


DEBIT ATION OF RECIPROCITY RELATIONS 


According to reference 5, the reciprocity relation for 
general three-dimensional unsteady motion can be written 

JJJ ~ (xi,yi,t^Wi{x u y u U) dxidyidti 

I* 

= JJJ — (x2,y2,h)W 1 (x 2l y2,t i ) dxdyifai (Cl) 

where the volinne of integration V is that swept out in x,y,t 
space by the wing. The subscript 1 refers to the wing moving 
in the forward direction and subscript 2 refers to the wing 
moving in the opposite direction in the same manner. The 
coordinate systems are related by 

X!=— x^+c— MT 


2/i= — y 2 +2s 

ti— — U.-\-T 

where s, c are wing semispan and chord, respectively, and T 
is some fixed value of time. These quantities are elucidated 
in figure 18. 


Substitution of these results into equation (Cl) yields 

Equation (C2) can be differentiated with respect to T, 
yielding 

The binomial expansion is now performed : 

(C3) 



Now lot the wing associated with the subscript 1 have the 
vertical velocity distribution 

’^K aJ r a )‘( £ ?)' 


In equation (C3) the spanwise integration is carried over 
the whole wing, but it can easily he reduced to integration 
over, say, the left panel by use of the factor [1 ( — l)* +rt ]/2. 

Thus, equation (C3) can he written 


By comparison with equations (36) and (37), it is seen that 
the integral terms in the last equation correspond to the 
generalized indicial forces /jjjjand //£, so that the summations 
can be written 


£<-i)"0/S-£<-iy 0/8 (04) 

where the quantity (g+n) must he an even number. 
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TABLE I.— VALUES OF GENERALIZED INDICIAL FORCES, F* 

The generalized indicial force coefficient F% is defined by 
equation (36). It is the response for a mode shape having 
a unit amplitude 



and a loading induced by a unit value of vj/ U 0 , 



The table gives values of F)* against time (actually chord 
lengths traveled) for 

1=0 

j= 0,1,2 

72 = 0,1,2,3,4,5 
27=0,1,2,3,4,5 
M=l.l, 1.2 
A = 4 


TABLE I.— VALUES OF GENERALIZED INDICIAL FORCES, F% 
(a) 1=0; j=0; Hf =1.1 


3.670 
3.602 
3.660 
3.276 
3.200 
a 363 
3.317 
3.310 
3.22" 

a 662 
aB32 
1. 230 
2.627 
2.010 
6.267 
6.800 
0.328 


7.273 
7.273 
7.273 
7.273 
7.272 
7.271 
7.209 
7.312 
7.671 
1.101 
1.213 
X 20 
L03 
2 09 
X01 
L 77 
110 


1 

2 

3 

4 

5 

\ n 

•\ 

0 

1 

2 

3 

D 

5 


0 

n 

2 

3 

D 

5 


4.848 

7.273 

1L04 

19.39 


1630 

1848 

7.273 

11.64 

19.39 

3125 


4.848 

7.273 

11. 64 

19.39 

3125 

5118 

1047 

2.863 

7.280 

minim 

19.48 


KSLJ 

1848 

7.277 

1L 65 

19. 45 

33.39 


4.848 

7.273 

11.64 
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3133 

58.41 

2687 

4.805 

7.321 

11.77 

19. 73 
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4.848 

7.289 

1L 70 

10.50 
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4.847 
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12.13 

20.89 
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W.S7 

3.620 

6.037 
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■liil 
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2L78 
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Kl:7ol 

7.273 

1L90 
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37.05 

mm 

a 062 

6.0S0 
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13.70 

2100 
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41.62 


4.836 
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20.91 

3110 

72. in 

3.670 

5. 174 

a 251 

14. 24 

20.00 
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12.95 
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42 76 
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7.320 
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39.48 
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a 871 

5. 604 

■taisa 
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30.03 


1763 

6. 128 

1152 

1118 

26.29 

51. 17 


5.0P2 

7.685 

12.92 

2132 

44.38 

Jtfl, 05 

1. 105 

0.023 

19. 10 

18.48 

sc. oe 


1 942 
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«. 774 

16.62 

2178 

66.86 



1132 
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7.269 

1280 
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52. CO 

■ 

1424 

6.217 

10.45 

19. (n 

39.79 

85.03 


Efftl 

1301 

1143 

31.87 

6129 

145.8 

5.325 

8. 561 

16.92 

3280 

73.29 


1872 

0. 992 

iais 

22 96 

61.30 

117.9 


1738 

1145 

19.07 

3168 

84.89 

1G14 

6. 7f$ 

o. »ro 

lasi 

39.81 

92 06 


5.203 

7.694 

la 69 

27.63 

Cl. 00 

147.6 


7.225 

11.34 

21.20 

44.43 

10L3 

2419 

a i88 

1L 21 

22 03 

5L27 

12d.O 


6.659 

■:gVTil 

16.73 

3142 

7134 

1911 


7.895 

12.63 

24.44 

53.42 

127.9 

329.2 

7.631 

la 86 

20.79 

71. 74 

188.0 


6.338 

9.866 

19.23 

43.28 

108.0 

291.4 



14.67 

29.72 

6168 

174.9 

4S0.0 

8.372 

15. 76 

34.88 

86.41 

2320 


6.834 


21.75 

60.36 

129.2 

3522 



1115 

3162 

70.62 

2010 

5S11 

■M&lll 

17.38 

38-74 

■IkiSal 

200.3 


■Sgill 

1L87 


65. S4 

144.1 

400.7 


exqi 

IT. 06 

3190 

82 21 

232.4 

6C7.7 

11. 61 

10.39 


58 18 

1024 


11. 6-4 

19.39 


58.18 

1014 

186.2 

5 

1139 

3125 

58 18 

1014 

1812 

3315 

11. 0J 

19.20 


53.32 

1018 


1L 64 

19.39 


5125 

1017 

180.9 


19.39 

3125 

6120 

1016 

1S1G 

339.8 

11.62 

19.43 


58.72 
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11.64 

19.39 


5145 

104.4 
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3125 

68.27 
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3413 
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33.83 
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59. 16 
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62.10 
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11.64 
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6182 

1018 
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1L 64 

10.79 
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61 52 

12LS 


11.04 

19.39 
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2112 


1139 

3125 

5123 

1018 
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1 L 61 

19.90 
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06.37 
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1L 04 

10.39 
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19.39 

3125 

5153 

1114 

2117 
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11. 71 
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37.83 

6S. 58 

1313 
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3145 

6151 

12L3 
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1152 
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6121 
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217.4 

431.4 

1221 

2L 23 
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76.86 
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20.49 


69.02 
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20.49 

35.12 

6109 
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242.8 

4913 

iaoi 
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179.0 


12.99 
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319.8 


21.67 

37.16 

6151 
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14. 67 
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3414 
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27.81 

53.44 

112.4 
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10S6. 

IS. 11 

32.80 

74.62 
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59.25 

1215 
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29.83 

5L6S 

1011 

2210 

Wttmm 

1331 
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KiXiM 

216.9 

5617 


19.56 

3166 

67.98 

1516 

377.1 

990.2 


32.83 

67.49 

1111 

269.6 

667. 5 

1705. 

2239 

48.67 

114.3 

296.0 

818.9 


22.12 

3192 

82.24 

195.9 

51L0 

1430. 


37.17 

6165 

142.7 

3412 


254L 

26.72 


132.2 

351.0 

100L 


23.98 

42.81 

92.51 

226. 3 

606.9 

1746. 


4134 

7130 

1614 

3810 

muwm 

3097. 

28.13 

*7). IN 

146 4 

390 9 

1119. 

■ 

20 10 

46.78 



672.8 

1951. 


44.02 

8110 

1713 

4314 

nn. 

3461 
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TABLE I. VALUES OF GENERALIZED INDICIAL FORCES, Fj* — Continued 

(d) Z=l;i=l; M=l.l 


ThU 

e 

\ n 

0 

1 

2 

3 

4 

6 

IX 

0 

1 

2 

3 

4 

6 

l>\ 

0 

1 

2 

3 

4 

5 

0 

.055 

.11 

.22 

.33 

.44 

.670 

.780 

LO 

1.671 

22 

2.76 

3.C07 

6.6 

7.333 

11.0 

0 

]| 
1 T*< 

IS 

ii 

ii 

ii 

1. 212 
1.212 
1. 213 
1.212 
1.232 
1.268 
1.289 
1.313 
L415 
1.620 
L770 
2000 
2.169 
2101 
2721 
2897 
2990 


2 424 
24 29 
2442 
2490 
2663 
2 675 
2779 
2872 
3.201 
3.549 
4.461 
6.391 
6.111 
7.193 
8.784 
9. 738 
10.06 

3.879 
3.890 
3. 931 
4. (Ml 
1.228 
1.161 
4.681 
4.8S1 
6.662 
6.299 
8.287 
10.49 
12.23 
14.96 
19.09 
21.81 
22.29 

6.465 
d 493 
a 675 

as 66 
7.290 
7.821 
8.302 
8.761 
10.24 
11.88 
ia63 
21.97 
2a 48 
33.77 
46.2 7 
53.68 
63.94 


li 

E 

igrn 


Z 424 
2.426 
a 431 
a 453 
a 496 
a 564 
a 638 

a 701 

a 942 

a 192 
asm 
4.409 
4.855 
a 492 
a 395 
a 905 
7.144 

3.879 

3.886 

3.902 

3.970 

4.0S0 

4.235 

4.390 

4.637 

6.017 

5.634 

6.861 

8.225 

9.284 

ia8i 

13.07 

14.41 

14.88 

a 406 
a 482 
a 632 

a 719 

a 994 
7.358 
7.700 
a 013 

9.080 
10.23 
ia29 
ia63 
19.25 
2a 32 
2a 45 
3a 39 
3105 

1108 

3113 

1128 

1173 

12.41 

13.27 
1105 
1178 
17.20 
ia85 

27.27 
35.83 
4a 86 
5114 
7175 
8130 
85.11 




a 879 

a 8si 
a 888 
a 021 
a 982 
1089 
1203 
1299 
1674 
a 063 
a 021 
a 939 
7.625 
a 599 
9.977 
10.75 
11.12 

a 465 
a 474 

a 500 

a 604 
a 776 
7.022 
7.270 
7.488 
a 276 
9. 106 
1L23 
ia39 
15.03 
17.46 
2100 
2a 09 
2a 84 

1108 
1111 
1119 
1149 
11 94 
12.54 

iaio 
iaei 
15.38 
17.27 
2a 29 
27.71 
3196 
3a 51 
4a 34 
5160 

sa oi 

ia39 
19.47 
19.69 
2a 47 
2L62 
23.07 
2t®8 
25L 61 
2a 70 
34. 17 
4a 69 

ea82 

72.41 
9a95 
119.8 
110.0 
lit 7 

r-NpUpfOMH o 

3 

2.424 
2.42 1 
2.426 
2435 
2460 
2509 
2567 
2 016 
2811 

3.012 
3.484 
2 000 
4.209 

4.012 
6.160 
6.415 
6.620 

|| 

E 

K 

H 

1L0S 
1L 10 
11.14 
It 31 
1L 69 
1200 
1241 
1278 
14. 10 
15.49 
10.02 
22 62 
25.34 
29.39 
35.21 
38.04 
39.91 

19.39 

19.44 
1167 
20.07 
20.83 
21.86 
22.81 
23.68 
2a 69 
29.91 

38.44 
47.61 
64.77 
66.76 
82.24 
92:6(1 
95.12 

34.48 
34.61 
34.98 
3a 33 
38.31 
40.81 
43.10 
45. 22 
62.30 
00.03 
8L44 
106.9 
125.6 

157.2 

2oai 

240.2 
243.4 

4 

a 879 
a 879 

assi 

a 896 
a936 
4.015 
4.109 
4. 186 
4.601 
4.825 
6.589 
a 276 
am 
7.431 
a 312 
& 730 
9. 089 

a 465 
a 465 
a 468 
a 493 

a 561 

a 692 
a 849 
a 976 
7.604 
a 046 
9.335 
10.51 
1L 37 
12.54 
14.15 
iao4 
1^64 

11.08 
11.09 
11. 12 
11.19 

11.38 
11.65 
It 97 
12.23 
13. 2S 

14.38 

17.01 
19.64 
21.43 
24.30 
27.87 
2a 99 

31.01 

19.39 

19.42 

19.49 

19.78 

2a 28 

2a 96 
21.67 
2a 29 
24. 57 
2a 96 
aa 04 
39.18 
4a 85 
5a 72 
6a 69 
6d 64 
ea 73 

3148 
3166 
3179 
35.65 
3a 97 
3a 73 
40.41 
4193 
47.18 

6a 68 

67.03 

saw 

95.92 
1119 
14a 3 
161.2 
1611 

V. X 

m 

fk * 


a 464 
a 465 
a 468 
a 493 

a 561 

a 692 
a 849 
a 976 
7.504 
a 044 
9.326 
10.48 
1132 
ia43 
ia92 
1160 
1&24 

1L0S 
It 08 
IL09 
11.13 
11.26 
It 47 
It 74 
11.90 
12:80 
13.79 
16.00 
18.01 
19.40 
2L49 
24.25 
25.77 
2a 63 

10.39 

19.40 
19.43 
ia53 
19.86 
2a 37 

2a 02 

21. 38 

2a 20 

25.08 
29.63 
3108 
37.35 
4196 
4a 47 
sa 14 
ca&2 

3148 
3152 
34.64 
35.14 
35.97 
37.20 
sa 45 
39.54 
4a 55 
47.76 
58.42 
69.19 
77.34 
sa 34 

ioa7 
nag 
12a 8 

62.06 
62.20 
62.61 
61. 11 
6a35 

eass 

72.63 
76.23 
84.68 
01.52 
120. 8 
143.9 
170.7 
201.2 
25tl 
283.6 
293.3 

lias 

113.2 
1114 
liae 
1218 
13a 7 
140.0 
14a 7 
JG9.0 
19a 5 
26a 7 

.ms 

398.3 

495.9 
64a 7 
75a 7 

761.9 


TABLE I.— VALUES OF GENERALIZED INDICIAL FORCES, Ffc— Continued 

(o) Z=0;j=0; Af=1.2 


Ucl' 

0 

\ » 

0 

i 

2 

3 

4 

6 

IBS 

0 

1 

2 

3 

4 

5 

P 

El 

■1 

mm 

3 

ra 


0 

.00 
.12 
.24 
.36 
.48 
.616 
.0 
.8 
L0 
L 6 
2.0 
24 

3.0 

4.0 

ao 

0 

3.333 

3.293 

3.266 

3.187 

3.123 

aoso 
a 068 
a 061 
a 147 

3.282 

3.652 

asas 

a983 

4.221 

4.633 

4.894 

3.333 
a334 
a334 
a337 
a340 
a 3i6 
a 347 
a369 
3.621 
a706 
4.163 
4.661 
4.833 
6200 
6 681 
6 173 

1444 
1448 
1460 
1499 
1556 
1019 
4.054 
1710 
a 013 
5.355 
a 202 
a 975 
7.531 

a 260 
9.196 
iaoe 

a667 
a 079 

a 713 

0.836 
7.000 
7.211 
7.325 
7.460 
8.107 
8.821 
10.62 
12.31 
la65 
16. 19 
17.28 
18.99 

10.67 
10.70 
10.79 
It 12 
It 69 
1214 
1245 
12 77 
14.23 
16.83 
iaoi 
2a 94 
2a 92 

3a 91 

35.95 

39.72 

17.78 
17.86 
iao9 
18. 91 
20.09 
2L50 
22.31 

2a 10 

2a 49 
30.22 
40.08 
50.01 

57.60 
67.85 
80.76 

89.60 


a 333 
3.333 
a 330 
3.323 
a 312 
a 299 
a 292 
a 307 
a 432 

a583 

a 952 
4.274 
4.498 
4.788 

a 172 
aoo2 

4.444 

4.444 

4.444 

4.446 

4.446 
4.440 
4.446 
4.474 
4.672 
4.907 
a 484 
a 997 
a 360 
a 833 
7.453 
8.094 

a 667 

a 67i 
ass2 
6.723 
0 781 
0847 
6885 
6 956 
7.366 
7.834 
6993 
10.05 
10.81 
1L80 

iao8 

14.28 

ia67 
10.68 
10.73 
1689 
It 12 
11.39 
It 64 
It 73 
12 67 
12 72 
16 33 
16 78 

2a £6 

22 92 
25.93 
28.46 

17.78 

17.83 

17.96 
18.45 
19.18 

19.97 
20.44 
20.93 
22 17 

25.83 
3L 92 

37.97 
42 48 
48.48 
6007 
61.87 

30.48 
3a 60 
3a 97 
32 28 
34.17 
36 41 
37.69 
38.91 
44.39 
60.38 
66 04 
81.67 
93.50 

109.6 

120.7 

14a 8 


4.444 
4.444 
4.444 
4.443 
4.441 
4.437 
4.435 
4.469 
4.046 
4.807 
a 401 
a 867 
a 192 
a 612 

7.162 

7.771 

6667 
6667 
6667 
6067 
6067 
6667 
6667 
6707 
7.002 
7.352 
8.311 
8.974 
9.613 
10.21 
It 13 
1209 

10.67 
10.67 
10.69 
10.74 
10.82 
10. 91 

iasc 

1L06 
1L69 
12.41 
14. 19 
iasi 
ia97 
ia49 
20.40 
22.31 

17.78 

17.80 

17.87 
18.11 
ia40 

18.87 
19.09 
19.39 
2a 87 
22.53 
2a 63 
30.55 
3a 37 
37.08 
4L83 
4a 88 

3648 
30.65 
3677 
31.55 
3267 
3898 
34.72 
36.62 
39.18 
4820 
6846 
68 29 
76 68 
80.29 
9259 
1021 


0 

.06 

.12 

.24 

,36 

,48 

.MG 

.6 

.8 

L0 

Lfl 

ao 

a4 

ao 

10 

ao 

3 

a 667 
a 667 
a 667 
a 607 

a 666 
a065 

6.604 
a 703 
a 993 
7.335 
a 165 
8.892 
0.390 
10.05 
10.92 
1L85 

10.67 
iae7 
iae7 
10.67 
10.67 
10.67 
iao7 
ia73 
1L 20 
1L 76 
iai3 
1135 
16.21 

16.32 
17.79 

10.32 

17.78 

17.79 
17.81 
17.89 
iaoi 

18. 14 
ia22 

18.39 

19.40 
2a 5S 
2a 48 
2a ii 
2a oo 

30.47 
3a 67 
3a 71 

30.48 
30.61 
3a 62 
3L01 
3L67 
32.22 
32 53 
3a 06 
36 63 
38.29 
46.19 
6L61 

6a 28 

62 43 
70.31 
77.07 

5a 33 
5a 40 

5a 82 

65.13 
57.00 
59.18 
ea 41 
6L 76 
07.95 
74.77 
92.12 
ioa 7 
121.0 
137.4 

15a 1 

174.3 

94.81 
95.18 
9a 22 
99.96 
105.4 
11L 7 
116.3 
na9 

134.7 
16L9 

iga 9 

24L4 

27a 0 

320.2 

377.2 

417.8 

1 

10.67 
10.67 
10.67 
10.07 
la 67 
10.67 
10.67 
10.73 
1L 20 
1L75 

13. 10 
14.29 

la 12 

la 19 
17.60 

19. 10 

17.78 
17.78 
17.78 
17.78 
17.78 
17.78 
17.78 
17.89 
13.67 
19.60 
2L8S 
2a 91 
2a 34 
27.20 
29.63 
32.18 

3a 48 
3a 49 
3a 53 
30.66 
3a 84 
3L08 
31.18 
3t 46 
33 18 
35.16 
40.07 
44.62 
47.70 
6L87 
67.27 
62 43 

5a 33 
sa4o 
53.58 
54.23 
5a 10 
5a 23 
5a 84 
57.66 
6L89 
6d 62 
7a 44 
89.43 
97.42 
107.9 
12L4 

13a 0 

94.81 
9a 03 

9a 66 

97.89 

10L1 

104.8 

106.9 

109.3 
120.0 
13L9 
162.0 
190.8 
212.1 
24a 3 
27a 2 

304.3 

170.7 
17L3 
17a 1 
176 7 

189.1 

200.2 
206.5 
212 9 
2465 
2768 

349.8 
427.7 
486 5 

565.0 

665.1 

7365 

1 

17.78 

17.78 

17.78 

17.78 

17.78 

37.78 

17.78 
17.88 
18.67 
19.60 
2LS6 

2a 88 
2a 26 

27.07 
29.45 

31.07 

30.48 
30.48 
30.48 
3648 
3648 
S64S 
3648 
36 66 
3201 
sa 60 

37.51 
40.98 
4843 
46 61 
66 79 
55.15 

53.33 
53.35 
53.42 
eae3 
sa 94 
54.30 
54.50 
54.09 
57.97 
61.41 
69. 91 
77.61 
8a 13 
9a 34 
09.70 
108.7 

| 

170.7 

171.1 

172.1 
17a 1 
18L6 
I88.'2 
191.9 
isa 0 
2iao 
23a 0 

289.2 
3ia 0 
377.5 

427.4 
49a 7 

540.5 

sias 

311.4 

314.7 
32a 3 

34a 1 

302.8 

374.0 
38a 3 
434.7 

m 8 

G29. 6 
768.3 

872.0 

ioia 

im 

131a 
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TABLE I— VALUES OF GENERALIZED INDICIAL FORCES, F % — Continued 


(f) 1=0; j=l; Af=1.2 


w 

c 

\« 
l \ 

0 


2 

B 

B 

6 

\ n 

ff \ 

0 

1 

2 

3 

4 

5 

ff \ 

0 

1 

2 

3 

B 

■ 

0 

■9 

Hi 


2.222 




1 

1.667 

2.222 

B!K!il 

■HI 

mm 

ia24 

2 



a 333 

8.889 

ia24 

2a 07 

.00 


■n?: 

■tffi 

2.222 

■nctti 

a 343 

acic 


1.664 

2.219 




15.28 


2.219 

3.329 

a 329 


ia26 

2a 74 

.12 


BTrtr 

1.659 

2.219 


a 373 



L 657 

2.211 

a 325 

a 340 

a 944 

ia43 


2.211 

a 317 

a 318 

8.894 

ia32 

2a 08 

.24 



L 635 

2.209 

a 363 

a4&i 

9.357 


L 027 

2. 178 

a 299 

a 355 

a 097 

iafi7 


2.177 

a 267 

5.270 

^KYriT 

ia65 

27.86 

.36 


■TTw! 

L 590 

2.189 

bh uNii 

a 645 

9.881 


L679 


a 252 

6.367 

a 315 

ia77 


2.120 

3.183 

a 185 


ia87 

29.10 

.48 


■TTn 


2.153 

■TTfitl 

a 821 

■ 0.1.91 


1.612 

2.045 

a 179 

a 357 

9.542 

17.71 


2.040 


a 057 


ia20 

3a 07 

.645 


■Tltr 


2.124 

"■if-Titl 

^■TCTiT 

■nTsI 


L469 

L 994 

a 126 


9.6M 

ia2i 


L 987 

2.989 

4.966 

a 797 

ia33 

3L40 

.6 


Bn~ rr 

L 489 

2.123 

Bil !k?1 


1L 19 


1.450 

L 974 

a no 

a 366, 

a 810 

ia78 


tm 

2.959 

4. 944 

8.838 

iaco 

32.39 

.8 



mwsm 

2.214 


a 732 

BETul 


L 451 

1.997 

a 224 

a 727 

B1'£U 

2L65 


1. 978 

2.992 

a 099 

a 387 

ia28 

37.15 

1.0 


1.313 

L56S 

2.357 

hi ittii 

7.600 

ia20 



2.069 

a 409 

a 219 

12.17 

2a 09 


2.041 

BtliiSi 

a 375 

iai6 

20.30 

42.80 

1.6 



L 760 

2.778 

BTTTPl 

iao7 

2L46 



2.307 

a 967 

7.367 

ia88 

34.93 


2.248 

KfEll 

0.217 

12.37 

2a 38 

59.10 

2,0 


Bwm 

Bl 11 1 

3.220 

■iTTff-1 

12.68 

2a 28 


L 771 

2. 564 

4.557 

a 128 

ia8i 

4a 58 


2.469 

a 831 

7.110 

BCiriH 

32.71 

7a 74 

2.4 


HTfplTl 


3.66S 

Bil vti 

14.78 

33.79 



2.769 

a 022 

iaso 


54.16 


2.644 

4.134 

7.818 

■iU 

37.74 

90.80 

3.0 


BlCtrl 

BtViffil 

4.002 

■n iffii 

17.74 

4L 64 


■ISI 


a 686 

11.96 

27.36 

oa32 


2.896 

4. 671 

a 827 

ia is 

44.84 

110.9 

■to 


L 949 

BtVivrl 

4.773 

Bil Ttl 

21.80 

52.24 


2-325 

a 509 

a 647 

14.30 

3a 43 

82.82 


a 284 

a 230 

ia29 

22.82 

04. 05 

13a l 

0.0 


2.268 


a 518 

|jj Ifl 

2a 14 

60.14 


2.694 

4.063 

7.688 

iac2 

38.55 

9a 38 



mg] 

11.01 


03.05 

160.1 

0 

3 


5.333 

&8S9 

ia24 

2a 67 

47.41 

4 

■SSI 

a 880 

ia24 

2a 67 

47.41 

sa 33 

5 

8. 889 

ia24 

2a 67 

47.41 

86.33 

16a 2 

.00 


B.l ;V>; 

5.327 

8.882 

ia24 

2a TO 

47.64 


a 327 

a878 

ia23 

2a 67 

47.46 

8a 66 



ia22 

2a 64 

47.40 

8a 42 

165.6 

.12 


Bj) :)lr 

6.307 

8.861 

ia 24 

2a 80 

47.03 


mwjri 

a 845 

iai9 

2a 67 

47.63 

8a 23 



la 16 

2a 68 

47.40 


i5a 7 

.24 



5.227 

8.775 

15.24 

27.16 

49.41 


5.227 

a 711 

iaw 

2a 65 

4a 22 

8a78 


■:%j ( | 

14.93 

2a 30 

47.31 

Kma 

101,2 

.30 


Bn rfS 

5.093 

0626 

ia2i 

27.07 

6L 61 


a 093 

a 489 

14.77 

2a 56 

4a 04 

92.67 



14.66 

26.83 

47. 15 

saoo 

107.8 

.48 


■n 3 

4.907 

8.403 

mo 

2a 16 

54.12 


■ Im'J 

a 178 

14.38 

2a 34 

4a 83 

9a 88 


a 178 

14.02 

2a 13 

4a 70 

89.39 

na 4 

.645 


■1 ;-:r 

4.782 

a 246 

14.99 

2a 36 

5a 42 


JE w?:*! 

7.971 

14.10 

2a ii 


99.09 



iaeo 

24.05 

4a 20 

8a 84 

179.2 

.0 



4.734 

8. 205 

law 

2a 79 

5a99 


m 

7.889 

14.03 

2a 20 




KIM 

ia&2 

24.51 

4a 42 

muLim 

181.0 

.8 


Lit 

4.786 

a 446 

15.93 

31.69 

6a 12 



7.976 

14.43 

27.69 

5a 66 



7.074 

ia67 

2a 18 

4a 99 

99.61 

209.2 

1.0 


■j] Iff; 

4.955 

8.8S9 

17.19 

3a 13 

74.91 


4.948 

a 257 

ia 17 

2a 83 

61.77 

13a i 


a 253 

14. 16 

2a 46 

52.71 

110.3 

239.7 

1.5 


3.415 

5.516 

10.25 

2a 84 

4a 20 

>0>l:B 


a 492 

9.189 

17.45 

3a 05 

7a 16 

181.9 


ai72 

ia76 

3a 39 

63.55 


32a4 

2.0 


3.707 

a 122 

11.69 

24.69 

5a 84 

132.8 


a 064 

■ nJuB 

mss 

42.61 

97.54 

234.3 


10. 16 

17.48 

34.59 

74.98 

173.3 

419.5 

2.4 



aoo3 

12.84 

27.72 

64.28 

isa 9 


■ IVNII 

SjtjitB 

2L80 

47.77 

112.1 

27a 3 


10.93 

ia86 

37.91 

83.00 

mo 

404.0 

3.0 



7.299 

14.47 

32.02 

7a 20 

191.0 


7.187 

12.16 

24.65 

5a 10 

132.7 

33a 9 


12.05 

20.83 

42.05 

oa77 

mi 

699.8 

4.0 


■JjtvJ 

8. 348 

ia85 

3a 10 

92.70 

237.4 


a 186 

EETjB 

2a 56 

65.48 

101.2 

E3XB 


ia74 

2a 82 

49.69 

114.9 

285.3 

74a 8 

ao 


im 

9.667 

1^49 

44.02 

107.0 

273.6 


9.483 

iaio 

3a 04 

7a 66 

isa 0 

4Sa4 


ia 9i 

27.68 

67.37 

132.8 


867.0 
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TABLE I.— VALUES OP GENERALIZED INDICIAL FORCES, Fg— Concluded 

(h) 1=1; j—l; M=1.2 



















































